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Nonlinear I-V characteristics have been observed in in- 
sulating quench-condensed films which are locally supercon- 
ducting. We suggest an interpretation in terms of the en- 
hancement of conduction by the depinning of a Cooper pair 
charge density wave, Cooper pair crystal, or Cooper pair glass 
that may characterize the insulating regime of locally super- 
conducting films. We propose that this is a more likely de- 
scription than the Coulomb blockade or charge-anticharge un- 
binding phenomena. 



In the context of the Bose-Hubbard model, which 
is equivalent to the model of a Josephson junction 
array with charging, zero-temperature superconductor- 
insulator(SI) transitions in two dimensions (2D), tuned 
by disorder or magnetic field, are believed to be direct, 
with metallic behavior only at the quantum critical point 
101 . Recent experiments have suggested the existence of 
a significant metallic phase between the superconductor 
and insulator. The Stanford group reported this for nom- 
inally superconducting MoGe films in moderate magnetic 
fields at low temperatures, and conjectured that it was 
due to dissipation They later found "true" super- 
conductivity in low fields g. Long ago, metallic be- 
havior was reported in quench-condensed granular films 
over a range of thicknesses intermediate between those 
for which films were insulating and superconducting [Q, 
and it was found more recently in Josephson junction ar- 
rays 1^. Re-examination of the theory has involved the 
inclusion of aspects of percolation elaboration of the 
Bose-Hubbard model 0,^, and consideration of dissipa- 
tive Bose systems 1^. Das and Doniach explained 
aspects of Ref. 4 by extending the Bose-Hubbard model 
to high filling and including nearest-neighbor as well as 
on-site Coulomb interactions. Their phase diagram con- 
tains the possibility of an intervening Bose metal phase, 
or a direct transition, depending upon the relative mag- 
nitudes of the various Coulomb and Josephson coupling 
energies. The Bose metal contains free vortices and an- 
tivortices, prevented from Bose condensing by dissipa- 
tion. The insulator is a condensate of vortices, or in other 
language, a Cooper pair charge density wave (CDW). 
Phillips and Dalidovich |^ , using the original form of the 
Bose-Hubbard model, also demonstrated that there was a 
Bose metallic phase. This followed from a subtle cancel- 
lation in the expression for the conductivity between an 



exponentially small population of bosonic quasiparticles 
and their associated exponentially long scattering time. 
The Bose insulator in this picture results from dissipative 
processes such as coupling to a heat bath, or from a high 
enough level of disorder. Ng and Lee using a different 
approach to theory of dissipative Bose systems, suggest 
that a Bose metal does not exist at zero temperature, but 
note a crossover regime at finite temperature that could 
be identified as a metallic phase. 

In this letter we shift the focus to the study of the in- 
sulating regime of granular quench-condensed films. In 
the past the insulator has not been carefully studied, as 
the issue considered was the condition for the onset of su- 
perconductivity. Since local superconductivity has been 
demonstrated in such films when they are insulating , 
they should be good candidates for the application of 
the Bose-Hubbard model, in which there is a nonvan- 
ishing order parameter amplitude. We have studied Ga 
films, which earlier exhibited the most striking metallic 
regime [|| . We propose that the nonlinear I-V character- 
istics found in the insulating regime may be evidence of 
a Bose insulating state, which could be either a Cooper- 
pair CDW or Cooper pair crystal. For sufficient disor- 
der it would be a Cooper pair glass. The nonlinearities 
appear to be similar to those associated with the depin- 
ning of charge density waves The behavior of these 
clustered films on the insulating side of the SI transition 
is different from that of nominally homogeneous films, 
where the observation of an orbital magnetoresistance 
linear in field was interpreted as evidence of vortices . 

Gallium films were grown on glazed alumina substrates 
held at liquid helium temperatures. This was done in an 
ultra high vacuum apparatus in which cyclic evaporations 
and in situ transport measurements at temperatures 
down to 150 mK could be made. Quoted thicknesses 
are nominal, being derived from a calibrated thickness 
monitor. Standard four-probe current-biased resistance 
measurements were made using a low noise dc current 
source and a nanovoltmeter. All electrical leads enter- 
ing the vacuum chamber were filtered with low-pass in- 
terference filters in series with 20-nH inductors. This 
combination exhibited voltage attenuations of 15 dB at 
1 MHz, increasing to 75 dB beyond 100 MHz. Magnetic 
fields could be applied up to 12 kC perpendicular and 20 
kG parallel to the plane of the sample. This apparatus 
is unique in that the evaporants are derived from com- 
mercial Knudsen cells, and the geometry of the cells and 
the source-to-substrate distance (60cm) is such that the 
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flux density at the substrate is uniform to better that one 
part in 10^. This permits the observation of the dramatic 
effect of very tiny changes in thickness, which might be 
averaged over in film growth systems with less uniformity 
of flux density. 

To set the stage, we show in Fig. 1 the evolution of 
R{T) with thickness for a series of a-Ga films, adapted 
from Ref. Q. The thinnest films, shown at the top, ex- 
hibit a monotonically increasing i?(T) with decreasing 
temperature. As thickness is increased, R drops near 
the bulk Tc of a-Ga, providing evidence for local super- 
conductivity. At lower temperatures it increases expo- 
nentially with decreasing T, until about 2K, where it 
levels off. Further increase in thickness results in metal- 
lic behavior at low temperatures. For somewhat thicker 
films R{T) ^ exp(^), which has been interpreted as 
evidence of quantum tunneling of vortices [^. At the 
lowest temperatures for these films, R(T) is also temper- 
ature independent. Further increase in thickness results 
in superconductivity. These changes with thickness are 
accompanied by changes in the /— ^characteristics at low 
temperatures. Metallic films have linear characteristics, 
whereas all of the other films have nonlinear ones. From 
the onset of local superconductivity to full superconduc- 
tivity, nominal film thickness increased from 12.75A to 
14.50A. 

The leveling off of R{T) of high resistance films was 
current dependent. The lower the current, the higher 
the value of R{T) that would be attained with decreas- 
ing temperature. This is shown in the temperature evo- 
lution of the I~V characteristics in Fig. 2. The differen- 
tial conductance, G{V), also shown in Fig. 2, exhibits a 
threshold voltage, Vr, which is independent of tempera- 
ture, and which decreases with increasing thickness, both 
of which are shown in Fig 3. 

Heating cannot account for the observed effects for sev- 
eral reasons. Considering the film to be a thermometer, 
one can define an "effective temperature," T^ff^ based 
on R{T). For a range of dissipated power spanning three 
orders of magnitude, from 10"^'^ to 10~^° W, T^ff is well 
above the thermometer temperature. However it did not 
increase with applied power, and in some cases even de- 
creased, implying that the changes in R(T) are not due 
to heating and that the film is not acting as a thermome- 
ter. Furthermore, the power at which the resistance lev- 
eled off as films became thicker decreased even when the 
fractional change in thickness was small. Actually, the 
power needed to heat to a particular temperature should 
increase rather than decrease with increasing film thick- 
ness, as both heat capacity and thermal conductance be- 
come larger. All of this implies that the nonlinearities 
are intrinsic. In addition, much higher resistances were 
measured in the same apparatus at lower temperatures 
for other films not exhibiting either superconductivity 
or flattening |p^ . Finally, the temperatures at which 
R(T) became temperature-independent were well above 
IK, making a scenario of not cooling the electrons un- 



likely Q. 

Quasiparticle tunneling involving superconductor- 
insulator-superconductor (SIS) junctions in the film 
[Qcan be ruled out because of the temperature indepen- 
dence of Vt shown in Fig. 3. If Vr were a supercon- 
ducting gap, it would increase with decreasing temper- 
ature. Not shown is the magnetic field independence of 
the threshold, which is also consistent with it not be- 
ing attributable to quasiparticle tunneling. However, the 
subthreshold conductivity does increase with increasing 
temperature and magnetic field suggesting the existence 
of quasiparticle transport channels. 

A phenomenological picture of charge unbinding has 
been used to explain nonlinear I — V characteristics in 
very much thicker granular Al films, driven insulating by 
a magnetic field In this picture, a charge-anticharge 
pair is formed on neighboring islands by tunneling, and 
the members of the pair attract each other. This attrac- 
tion can be overcome by an applied electric field which 
creates free charges. This model, which includes five free 
parameters that cannot be independently measured, can 
be fit to our data. The fit, which is not shown, does not 
reproduce the sharpness of the conductance increase at 
Vt- 

In Josephson junction arrays with ultra-small junctions 
|]l8| and Pd films (l9|, nonlinear behavior at high resis- 
tances has been attributed to effects related to the charge 
Kosterlitz Thouless Berezinskii (KTB) transition. As a 
consequence of the electrical duality between the charge 
and vortex KTB transitions [Q, nonlinear I-V charac- 
teristics would be expected, with / cx V"", and a exhibit- 
ing a jump from 1 to 3 at the transition temperature. 
There should also be a square root cusp in R{T) above 
the KTB transition temperature. These effects were not 
seen. This model is most likely irrelevant because the 
screening length in these films is too short for the inter- 
action to be logarithmic over a substantial distance, a 
requirement for this scenario pT[ |. 

Nonlinearities in the I-V characteristics could be the 
signature of a Coulomb blockade |^^. The fact that the 
high-voltage current asymptote of the I-V characteris- 
tic, at low temperatures, extrapolates back through the 
threshold voltage Vr, rather than through the origin, is 
consistent with this picture. The macroscopic I-V char- 
acteristics of a film would then be dominated by those 
associated with tunneling in and out of a single "repre- 
sentative island," according to the arguments of Ambe- 
gaokar, Langer and Halperin |2^]. The threshold voltage 
would be given by the energy necessary to overcome the 
Coulomb energy Eq in adding charge to this island, as 
Vt — Ec/e = e/2C , where C is the capacitance of the 
island. The central argument against this picture is that 
nonlinearities are not seen until the temperature corre- 
sponding to the onset of local superconductivity. Values 
of Vt ranging from 6mV to 130mV in films of differ- 
ent thicknesses, correspond to Coulomb energies much 
greater than ksT at 8.3K, where the threshold appears. 
If the nonlinearities were a Coulomb blockade effect, they 
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would appear at higher temperatures and would not re- 
quire local superconductivity as a precondition for their 
observation. 

We now turn to a picture of the insulating state in 
which it is CDW-like configuration of Cooper pairs or a 
two dimensional Cooper pair crystal. This idea is men- 
tioned in theoretical papers, but is not treated in detail 
10,0-^. We suggest that the nonlinearities are associ- 
ated with the depinning of this putative CDW, crystal 
or glass state. Features of the depinning that qualita- 
tively resemble the data are the sharply enhanced con- 
ductance at Vt, and the fact that at low temperatures, 
all the I-V curves come together at high bias with the 
asyrnptotic limit of the current extrapolating back to 
Vt Atomic force microscope (AFM) pictures of 

the film sequence discussed in detail here, which are not 
shown, were obtained at a thickness of 18.24 A, and af- 
ter the film was warmed up and removed from vacuum. 
Although on warming, annealing of the microstructurc 
takes place, the mesoscale structure should be largely un- 
changed, because room temperature, or 20°C, is not high 
enough to result in atomic diffusion over distances suf- 
ficient to change the mesoscale structure. Clusters form 
during growth, with the kinetic energy of the imping- 
ing atoms making them mobile on the substrate surface, 
even though it is attached to a holder cooled to helium 
temperature. The AFM pictures showed flat clusters of 
200 A diameter, which should be representative of the 
structure at low temperatures. Similar studies of a thin- 
ner film from a different series, revealed smaller clusters, 
suggesting a correlation between cluster size and thick- 
ness. It may be possible to consider the film structure 
of a random array of nearly identical clusters to be the 
origin of the pinning of the Cooper pair CDW or crystal. 
The decrease of the threshold voltage with film thickness, 
shown in the inset of Fig. 3, would then result from a 
reduction of the areal density of pinning centers with in- 
creasing thickness and physical cluster size. Scaling of 
the I — V characteristic near Vt , which is found in many 
CDW systems could not be tested. This was a conse- 
quence of our measurements being current biased. At Vt 
the voltage increased very rapidly, limiting the data in a 
manner that made a definitive test of scaling impossible. 
What is most important is that the nonlinear behavior 
manifests itself only when a local minimum in R(T) is 
present. Thus the phenomenon must be associated with 
Cooper pairing. 

The behavior in magnetic field should be noted. In 
fields up to 12 kG applied perpendicular and up to 20 
kG applied parallel to the film, Vt is unchanged, but the 
onset temperature for local superconductivity is lowered 
by up to IK. Field increases the conductance for V ^ Vt- 
Neither maximum value of field was high enough to com- 
pletely obliterate the threshold. The fact that Vt is field- 
independent suggests that the pinning mechanism is as- 
sociated with phenomena unaffected by field, i.e., defects. 
The effect of a parallel field on the subgap conductance 
is smaller than that of a perpendicular field, but is of the 



same order of magnitude. For a 17.2A thick Ga film at 
AK, a 12kG exhibited a magnetoconductance of -1-54% in 
a perpendicular field and -1-20% in a parallel field. The 
parallel field effect cannot be explained by misalignment 
with the plane, which was at most I*'. These observa- 
tions suggest that short range effects of magnetic field on 
individual grains are comparable to longer range effects 
probed by a perpendicular field. Alternatively the depen- 
dence on parallel magnetic field could imply some role for 
spin degrees of freedom. The sub-threshold conductance, 
as mentioned, most likely comes from quasiparticles. The 
application of a magnetic field would increase the number 
of such quasiparticles participating in transport. 

The picture we have presented here as well as the 
Coulomb blockade and charge-unbinding effects all in- 
volve Coulomb interactions. The difference between 
charge-structure depinning and these other phenomena is 
whether the behavior is collective, or whether charges in- 
dividually overcome a Coulomb energy. Further support 
for collective behavior is the similarity and even electrical 
duality of G{V) characteristics of "high" resistance films 
with R{I) characteristics of low resistance films. In the 
latter, the current threshold is a consequence of the de- 
pinning of vortices |Q. The saturation of G{V) and R{I) 
at low temperatures would imply that both types of de- 
pinning may be quantum phenomena. If collective charge 
behavior indeed governs the insulating phase of these 2D 
superconducting systems, then it is a member of a larger 
class of problems including vortex lattices, charge density 
waves, Wigner crystals, and magnetic bubbles, in which 
it is necessary to understand the dynamics of a periodic 
structure in the presence of static substrate disorder. It 
would then be appropriate to bring to the problem the 
ideas used to study such solids |^ . 
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FIG. 1. Evolution of the temperature dependence of the 
sheet resistance for a series of Ga films, adapted from Ref. 4. 
Eilm thicknesses range from 12.75 A to 26.67 A and mcrease 
from top to bottom. Thickness was changed in increments of 
O.O5A. Some data has been omitted for clarity. 



FIG. 2. (a) The evolution of current as a function of voltage 
with temperature in a 18.2 A Ga film. The temperatures are 
8, 7, 6, 5, 4, 3, 2, 1.5, 1, 0.5, 0.15 K from top to bottom, with 
curves overlappping at 2K and below, (b) Differential conduc- 
tance, dl/dV, as a function of voltage for the temperatures on 
the same film. A temperature independent threshold voltage 
of about 6 mV can be seen where the conductance suddenly 
increases. 



FIG. 3. Threshold voltage as a function of temperature of 
an 18.2 A Ga film. Inset: threshold voltage as a function of 
thickness for four Ga films in a series. 
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